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Abstract 

The neutron-rich carbon isotopes 19 > 17 C have been investigated via proton inelastic 
scattering on a liquid hydrogen target at 70 MeV/nucleon. The invariant mass 
method in inverse kinematics was employed to reconstruct the energy spectrum, in 
which fast neutrons and charged fragments were detected in coincidence using a 
neutron hodoscope and a dipole magnet system. A peak has been observed with an 
excitation energy of 1.46(10) MeV in 19 C, while three peaks with energies of 2.20(3), 
3.05(3), and 6.13(9) MeV have been observed in 17 C. Deduced cross sections are 
compared with microscopic DWBA calculations based on p-sd shell model wave 
functions and modern nucleon-nucleus optical potentials. J n assignments are made 
for the four observed states as well as the ground states of both nuclei. 
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With the advent of new radioactive beam facilities capable of producing in- 
tense beams of various nuclear species far from the stability line, even at 
the drip-line, an increasingly large amount of information on nuclear levels 
and modes of excitation is being accumulated throughout the nuclear chart. 
New phenomena, such as nuclear halos and skins PQ, enhanced El transition 
strengths [215] . and modifications of shell closures [U5j have been revealed. Of 
particular interest in recent years are the neutron-rich carbon isotopes, which 
have attracted attention not only from their own structural interest, such as 
the anomalously reduced E2 transition strength in 16 C |5f7] . but also for their 
implications for a new magic number at the neutron number N=16 proposed 
in oxygen isotopes [5118] . 

Email address: satou@phys.titech.ac.jp (Y. Satou). 
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Neutron number dependence of ground state deformations of carbon isotopes 
has been investigated in a deformed Hartree-Fock (HF) theory [9]. Generally 
the prolate deformation is expected at the beginning of the shell whereas oblate 
deformation arises towards the end of the shell. Note that the occurrence of 
nuclear deformations in the ground state is a consequence of the spontaneous 
symmetry breaking effect known in many fields of physics, and there is broad 
interest in elucidating the intriguing deformation-driving mechanism in atomic 
nuclei [10] , which would be a nuclear physics analogue of the Jahn- Teller effect 
in molecular physics The theory predicts prolate deformations for carbon 
isotopes with iV=9-ll. For 19 C with iV=13 two almost degenerate deformed 
minima are predicted with spins J n =l/2 + (prolate) and 3/2 + (oblate). Since 
the shape change at a neutron number smaller than the middle of N=8 and 20 
indicates a new shell closure at N=16, it is argued that definite information 
on the structure of 19 C is important to clarify the possible new shell effect in 
the neutron-rich carbon isotopes [9lfT2] . 

The 19 C nucleus is the heaviest odd carbon isotope. It is loosely bound with a 
neutron separation energy of >S n =0.58(9) MeV [T3], and exhibits one-neutron 
halo structure as evidenced by large Coulomb break-up cross sections p3] . The 
ground state spin and parity were also investigated via the measurements of 
longitudinal momentum distributions of charged fragments after the removal 
of one neutron from 19 C [ToTTE] . From these measurements the spin-parity was 
assigned to be J^ s =l/2 + . Other possibilities, however, have been suggested 
by authors in Ref. p2] based on their experiment in search of an isomeric 
transition. Moreover, there remains a controversy over the different widths of 
longitudinal momentum distributions measured at different energies [HUH], 
which has led to a conjecture of a possible resonance state just above the 
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particle decay threshold as a clue to solve such an inconsistency [20J. 

In this situation it it worthwhile accumulating experimental information on 
the ground as well as excited states of 19 C. This paper reports a new measure- 
ment using the (p,p') inelastic scattering on 19 C employing the invariant mass 
method in inverse kinematics. Both decaying neutrons and charged fragments 
were detected, and an isolated level was identified in the final state. The mea- 
surement was also made, partly for calibration purposes, on another loosely 
bound nucleus 17 C, having S n =0.73(2) MeV p3], for which J^ s is consistently 
reported to be 3/2+ [T6ll2Tll22ll23] . 




The (p,p') reaction has the following advantages: (1) The magnitude of the 
cross section, being sensitive to both the initial and final state wave functions, 
is configuration dependent. (2) The shape of the angular distribution also 
depends on the configuration through characteristic transferred L dependences 
of partial amplitudes. (3) Theoretical methods such as the distorted-wave Born 
approximation (DWBA) can be used to provide a first interpretation of data. 
There exists one (p,p f ) work on 19 C reporting two bound states at 0.20 and 
0.27 MeV using 7-ray spectroscopy [21]. These states are tentatively assigned 
as 3/2 + and 5/2 + , respectively, based on the assumption of J* s =l/2 + . Close 
proximity of levels near the ground state has made it difficult to identify levels 
from comparisons in excitation energy between theory and experiment. Note 
that in shell model calculations [16] the triplet of levels are predicted below 
0.62 MeV with no ground state configuration favoured. We have chosen to 
probe states in the unbound region, where in contrast to the bound region 
a lower level density is predicted up to about 3 MeV from the threshold in 
a shell model calculation described later. With no states known above the 
particle decay threshold, the measurement involved a search for resonances in 
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this region, and we report a new state in this paper. For C eleven new states 
up to 16.3 MeV excitation energy have been recently reported from the three- 
neutron transfer reaction 14 C( 12 C, 9 C) 17 C, with limited spin assignments |25j . 
This nucleus is also known to have a low- lying triplet of levels l/2 + , 3/2 + , 
and 5/2+ below 0.33 MeV [23] . 

The experiment was performed at the RIKEN Accelerator Research Facil- 
ity (RARF). The radioactive beams of 19 C and 17 C at 70 MeV/nucleon were 
produced using the projectile-fragment separator, RIPS [25], from a 22 Ne pri- 
mary beam at 110 MeV/nucleon. Typical beam intensities were 260 cps for 
19 C and 10.4 kcps for 17 C with momentum spreads AP/P of 3.0% and 0.1%, 
respectively. The secondary target was a cryogenic hydrogen target [27J hav- 
ing a cylindrical shape, 3 cm in diameter and 120±2 mg/cm 2 in thickness. 
The target was surrounded by forty-eight Nal(Tl) scintillators used to de- 
tect de-excitation 7-rays from the charged fragments. Each crystal had a di- 
mension of 6.1x6.1x12.2 cm 3 . The charged fragments were bent by a dipole 
magnet behind the target and were detected by a plastic counter hodoscope. 
Multi-wire drift chambers placed before and after the magnet were used to ex- 
tract trajectory information of the charged particles. Neutrons were detected 
by a neutron hodoscope consisting of two walls of a plastic scintillator ar- 
ray placed 4.6 and 5.8 m behind the target. Each wall had a dimension of 
2.14 w x0.72 H (or 0.90 H )x0.12 T m 3 . The total efficiency of the hodoscope was 
24.1±0.8% for a threshold setting of 4 MeVee. This was deduced by measuring 
the 7 Li(jo, n) 7 Be(g.s.+0.43 MeV) reaction at E p =70 MeV and using existing 
cross section data [28]. The invariant mass of the final system was calculated 
event-by-event from the momentum vectors of the charged particle and the 
neutron. A series of studies probing unbound resonance states in beryllium 
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isotopes has been successfully performed using a detector setup similar to the 
present one [29)30J31j . 

Figured] shows relative energy spectra for the (a) 1 H( 19 C, 18 C+n), (b) H( C, C+n), 
and (c) 1 H( 19 C, 16 C(2 + ; 1.77 MeV)+n) reactions, integrated over center-of- 
mass angles up to c . m .= 64°. The effect of the finite detector acceptance was 
corrected for in panels (a) and (b), i.e. in inelastic channels. Shown in the insets 
of panels (a) and (b) are spectra obtained in coincidence with de-excitation 
7-rays from the first 2+ states at 1.58(1) MeV in 18 C [32J and at 1.77(1) 
MeV in 16 C [23], respectively. The spectrum in panel (c) also required the 
coincidence detection of 7-rays from the 2 + state in 16 C. Background contri- 
butions from various window materials surrounding the target, measured with 
an empty target, are subtracted. Error bars represent statistical uncertainties. 
Peak structures at about 0.9 MeV in relative energy in 19 C and about 1.5 MeV 
in 17 C are clearly seen in spectra without requiring 7-ray coincidences. They 
are absent in the respective 7-ray coincidence spectra. We thus conclude that 
charged particles resulting from the decay of these peaks are in their ground 
states. In panel (b) we can see two more peaks at about 0.6 and 3.6 MeV 
in relative energy. These are visible in the 7-ray coincidence spectrum in the 
inset of panel (b); the 0.6 MeV peak is more clearly populated in the spectrum 
in panel (c). This observation indicates that these peaks in 17 C decay through 
the first 2+ state in 16 C. 

The experimental spectra were analyzed to extract the resonance energy E r 
and the width r r in the following way. Firstly, a single Breit-Wigner shape 
function was generated for certain values of E r and r r : 

T/ (.£/ re ])r r 

a{Elcl) ~ {E r + ae^) - E TcX y + {r^eO/2} 2 - (1) 
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The shift function AEi(E Te \) and the level width Ti(E re \), which depend on 
the relative energy E re \, were calculated by using the penetration Pi and shift 
Si factors [34] by the relations: 



AEi(E Tel )=T r x {Si(E r ) - S l (E Tel )}/{2P l (E r )}, 
Ti(E rcl )=T r xPi(E rcl )/Pi(E r ), 



(2) 
(3) 



where / refers to the decay angular momentum. Decay neutrons were supposed 
to be in the 1=2 orbit. The channel radius was taken to be R=r (Al/ 3 + A 1 / 3 ) 
with tq=1A fm, where A n and Af are the neutron and charged fragment mass 
numbers, respectively. Then, the experimental resolution, including the detec- 
tor resolution, beam profile, Coulomb multiple scattering of charged particles, 
and their range difference in the secondary target, was incorporated to gener- 
ate a response function. The relative energy resolution was well simulated by 
AE TC \=0. 13^ Evci MeV (rms). The procedure was repeated for each resonance 
peak by varying the initial values of E r and T r . Parameters which gave the 
best fit to the data were obtained by minimizing the x 2 . The solid curves in 
Fig. CD show the results of the fit; the dashed lines background introduced to 
reproduce the overall spectrum. 

Resonance parameters extracted are summarized in Table [TJ The excitation 
energy E x is obtained by the relation: E x =E r +S n +E* , where E* refers to the 
excitation energy of the daughter nucleus. The state at E x =l. 46(10) MeV in 
19 C was observed for the first time in this measurement. The 2.20(3), 3.05(3), 
and 6.13(9) MeV states in 17 C are respectively close to the 2.06, 3.10, and 6.20 
MeV states observed in the three-neutron transfer reaction [25]. Extracted 
widths for the 1.46 and 2.20 MeV states are T r =0.29(2) and 0.53(4) MeV, 
respectively. These values are of the order of the respective single-particle 
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estimates for the width of 1=2 resonances |35j, r sp =0.25 and 0.78 MeV, sup- 
porting the assumption of 1=2 decays of these states. For the 3.05 MeV state 
the T r value was not extracted since it was insensitive to the spectrum shape 
owing to the very small penetration factors at low relative energies. Cross 
sections leading to these states, angle-integrated up to 6> cm .=64°, are given 
in Table [TJ Differential cross sections leading to the 1.46 MeV state in 19 C 
are shown in Fig. [2, and those to the 2.20 and 3.05 MeV states in 17 C are 
shown in Fig. [3j The errors shown are the quadratic sum of statistical and 
systematic uncertainties; the latter in the absolute magnitude of the cross sec- 
tion is estimated to be 7%, including ambiguities in target thickness, neutron 
detection efficiency, and fitting procedure. The angular resolution varied from 
A0 c . m .=3.4° (3.1°) in rms at # c . m =4° to 5.1° (4.6°) at 60° for 19 C ( 17 C). 

To obtain an understanding of the results shell model calculations were per- 
formed using the code OXBASH [SB] within the Qftw configurations in the p-sd 
shell model space. For 19 C, the second 5/2^ state is predicted at 1.40, 1.54, and 
1.47 MeV with the WBT, WBP [37], and MK [38] interactions, respectively. 
These energies are close to the experimental value of 1.46 MeV; it is likely 
that this state corresponds to the 5/2^" state. The closest higher- lying state 
is the 7/2f state predicted 0.7-1.0 MeV above the 5/2^ state. For 17 C, with 
the WBT interaction four states 5/2^(1.72 MeV), 7/2^(2.33 MeV), 9/2^(3.01 
MeV), and 3/22~(3.08 MeV) are predicted above the decay threshold and below 
3.5 MeV. 

In order to further clarify the nature of states involved in transitions shown 
in Table [TJ microscopic DWBA calculations were performed using the code 
DW81 [39]. The optical potential was taken from the global parameteriza- 
tion KD02 (30] • A microscopic optical potential (UJ based on the approach of 
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Jeukenne, Lejeune, and Mahaux (JLM) [12] was also tested. The projectile- 
nucleon effective interaction was the M3Y interaction [43] . The transition den- 
sity was calculated with the shell model using the WBT interaction [37]. The 
single-particle wave function was generated in a harmonic oscillator well. The 
oscillator parameter was chosen so that the rms radius corresponding to the 
wave function reproduces the experimental value [44J: 6=2.07 fm for 19 C and 
1.83 fm for 17 C. The effect of core polarization on quadrupole transition ampli- 
tudes was taken into account in the isoscaler channel, by introducing isospin 
dependent polarization charges obtained in the FH+RPA particle-vibration 
model and parameterized in Ref. [10]: 5t=o=0.17 for 19 C and 0.22 for 17 C. In- 
tegrated DWBA cross sections are given in Table [1] for each transition listed. 
The J n values are discussed below. 

In Fig. [21 the DWBA predictions of the differential cross section leading to 
the 1.46 MeV state in 19 C are shown for three possible ground-state configu- 
rations. The solid line was obtained by using the KD02 optical potential for 
the supposed 1/2^^5/2^ transition, while the dot-dashed line by using the 
JLM potential for the same transition. The latter was specifically derived for 
the p+ 19 C system at E p =70 MeV. These curves agree closely with each other, 
justifying an extrapolated use in mass number of the KD02 potential (the 
nominal range is A=27-209) down to the A=19 region even for nuclei with 
large neutron/proton ratios. We therefore adopt the KD02 potential below. 
Dashed and dotted curves respectively assume J^ s =3/2+ and 5/2$, and the 
same final state. These curves fail to reproduce the magnitude of the cross 
section. Clearly, the data are much better described by the solid and dot- 
dashed curves obtained with JJ s =l/2f . If the 1.46 MeV levels is not 5/2$, 
it might be identified with the 7/2^ state. DWBA cross sections exciting this 
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7/2 + state from any of the members of the low-lying triplet, integrated up to 
6c.m.=64°, are, however, at most 4.2 mb (for 5/2^—^7/2^), only 50% of the 
data. We thus conclude that the spin of the ground state of 19 C is consistent 
with J^=l/2+ [MT5IT6] . and that of the 1.46 MeV state is 5/2+. 

In the study of the 14 C( 12 C, 9 C) 17 C reaction [25], a state was observed at 2.06 
MeV in 1 C, which was assigned as either 3/2 + or 7/2 + using the results of 
shell model calculations. Although there is a slight discrepancy in the exci- 
tation energies, a possible counterpart of this state is the 2.20 MeV state in 
the present (p,p') study, for which we examine the two proposed J 77 assign- 
ments. Of all possible transitions connecting the low-lying triplet of states 
and the (3/22", 7/2f) state, it is found that only three have sizeable cross 
sections when integrated up to # c m .=64°: 2.3 mb for l/2 + — >3/2 + , 2.7 mb for 
3/2+^7/2+, and 3.0 mb for 5/2+^7/2+. The DWBA cross sections for other 
transitions are at most 0.66 mb (3/2+— »3/2 + ). Of the three the l/2 + — >3/2 + 
case is readily excluded since it has been clearly demonstrated in a g-factor 
measurement that Jg s of 17 C is different from l/2 + [21]. We are thus left with 
the 7/2f assignment for the 2.20 MeV state. DWBA predictions leading to 
this state are shown in Fig. [3] for three initial-state configurations: solid line 
assumes J g w s =3/2 + , dashed one 5/2 + , and dotted one 1/2+ . The l/2f— »7/2f 
assumption gives much lower cross sections than the data. Of the remaining 
two we see that the 3/2+ assumption for the ground state better describes the 
data by reproducing the slope of the angular distribution, although the 5/2 + 
assumption also gives a reasonable description of the data. 

In Fig. [3j differential cross sections leading to the 3.05 MeV state are also 
compared with DWBA predictions obtained by assuming different configura- 
tions for the ground state, and the 9/2 + state for the excited state. This state 
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was populated only weakly. It had a sizeable cross section only at backward 
angles, where transition amplitudes with large angular momentum transfers 
are involved. It could be identified as the 9/2 + state reported at 3.10 MeV 
in the three-neutron transfer work [25], and predicted at 3.01 MeV in the 
present shell model calculations. The cross sections leading to this state are 
in good agreement with the solid curve calculated for the 3/2 + — >9/2 + transi- 
tion, excluding the possibility of J^ s =5/2 + for 17 C. This further corroborates 
the earlier 3/2 + assignments for the ground state of this nucleus ^611211122(23] . 
giving us confidence in the current procedure employing the shell model wave 
functions and DWBA calculations. 

The 6.13 MeV state in 17 C appears to correspond to the 6.20 MeV state 
reported close in energy in the three-neutron transfer work [23] . In that study 
it was tentatively assigned as either 5/2 J or 5/2 + , with a preference for the 
latter assignment due to larger occupancies predicted for the Od shells. The 
5/2 + and 5/2 + states are located at 6.25 and 6.81 MeV, respectively, in the 
present shell model calculations using the WBT interaction. DWBA cross 
sections leading to the 5/2 + and 5/2 + states, angle integrated up to 6 , c . m .=64°, 
are 0.94 and 0.25 mb, respectively. The experimental value of 0.8(1) mb favours 
the 5/2^" assignment for this state. The 5/2 + state is predicted to have a large 
occupancy of the 0d3/2 shell of 23%, which corresponds to almost one neutron 
in this shell, in contrast to other lower energy states involved in this study: 
4%, 4%, and 6% for the ground, 7/2+, and 9/2+ states, respectively. This 
indicates that promoting one neutron from lower-lying shells to the 0c?3/2 
shell is the dominant excitation process of this state. The deduced narrow 
width is consistent with the reported value of T r =0. 35(15) MeV in Ref. [25J. 

The 5/2 + shell model state at 1.72 MeV in 17 C is not observed clearly pre- 
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sumably due to the large expected width [25J. The angle integrated DWBA 
cross section leading to this state from the 3/2^ ground state was calculated 
to be 0.69 mb. 

According to the argument of Ref. [9] the spin of s =1/2^ for 19 C may imply 
a prolate intrinsic deformation for the ground state. In terms of the Nilsson 
diagram this indicates that four neutrons occupy the [220|] and [21l|] orbits, 
and one neutron occupies the [21l|] orbit in the ground state, where [An^Af^] 
are asymptotic quantum numbers referring to large prolate deformations. If 
one assumes that the observed l/2 + — >5/2 + transition is primarily associated 
with the promotion of the least bound neutron in the [21l|] orbit to the [202 1] 
orbit, one would obtain a value for the quadrupole deformation parameter of 
/5 2 ~ 0.4, by referring to neutron single-particle levels in a deformed Woods- 
Saxon potential |45j, and by equating the excitation energy of 1.46 MeV with 
the energy difference of the two Nilsson orbits. Interestingly, this value is 
consistent with the result of more recent deformed Skyrme HF calculations 
predicting a local prolate minimum with J n =l/2 + at /52=0.39 [10]. In the 
HF calculation, however, the ground state was predicted to be oblate with 
(32=— 0.36 and J 7r =3/2 + , and to be more bound by 2.05 MeV than the prolate 
minimum. Moreover a local oblate minimum having J 7r =l/2 + with (32=— 0.35 
was also predicted to be almost degenerate with the J 7r =3/2 + ground state. 
The present cross section for the l/2 + ^5/2 + transition observed in 19 C will 
be useful in distinguishing between the two l/2 + states with different signs of 
/?2, and provide a clue to further investigate the persistence of the new magic 
number A=16 proposed in oxygen isotopes [5|8] down to the carbon isotopes. 
The HF model [10] gives a good account of the ground state spin of 3/2 + for 
17 C. It is of interest to see if it also accounts for other states observed in this 
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experiment. 

In summary, we have demonstrated that the (p,p') reaction using the invari- 
ant mass method in inverse kinematics leading to unbound resonance states 
in the residual nucleus is feasible for structure studies even for nuclei far 
from stability. The measurements were made on 19 > 17 C at 70 MeV/nucleon. 
One resonance in 19 C and three in 17 C were observed above the particle de- 
cay threshold. A DWBA analysis employing shell model wave functions and 
modern nucleon-nucleus optical potentials was used to identify the transitions 
observed. The spin-parity of the ground state of 19 C was found to be consis- 
tent with l/2 + , and that of the strongly excited 1.46 MeV state was assigned 
to be 5/2t ■ F° r 17 C the observed states corresponded well with those reported 
in a three-neutron transfer study [25]. By adding information from this exper- 
iment spin-parity assignments of 7/2+ for the 2.20 MeV state and 5/24 for 
the 6.13 MeV state were made. The spectroscopic information from this study 
will impose stringent constraints on further theoretical investigations of light 
neutron-rich nuclei in this region. 

The authors acknowledge invaluable assistance of the staff of RARF during 
the experiment, particularly, Dr. Y. Yano, Dr. A. Goto, and Dr. M. Kase, and 
useful discussions with Professor I. Hamamoto and Professor J. A. Tostevin. 
This work was supported, in part, by the Grant-in- Aid for Scientific Research 
(Nos. 15540257 and 15740145) from MEXT Japan. 
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12 3 4 
Relative energy [MeV] 

Fig. 1. Relative energy spectra integrated over an angular range below 9 c , m =Q4° for 
the (a) 1 H( 19 C, 18 C+n), (b) ^("C.^C+n), and (c) 1 H( 19 C, 16 C(2+; 1.77 MeV)+n) 
reactions at 70 MeV/nucleon. Shown in the insets of panels (a) and (b) are spec- 
tra obtained in coincidence with de-excitation 7-rays from the first 2 + states of 
the respective charged fragments. The instrumental background was subtracted in 
the spectra. Solid lines represent the results of the fit, dashed lines background 
introduced to reproduce the overall spectrum. 
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Table 1 

Resonance parameters and populating cross sections from this experiment. The 
excitation energy and the cross section are compared to theoretical values. Cross 
sections are integrated up to 6> c . m .=64°. DWBA cross sections odwba leading to 
states specified by J w are obtained by using the shell model wave functions with 
the WBT interaction, assuming =1/2+ for 19 C and 3/2+ for 17 C. 



Nucleus E r 



Experiment 



'cxp 



Theory 



E x ffDWBA J n 



(MeV) (MeV) (MeV) (mb) (MeV) (mb) 



19 



C 0.88(1) 1.46(10) 0.29(2) 8.6(4) 



1.40 6.9 5/2^ 



17 



C 1.47(2) 2.20(3) 0.53(4) 3.8(2) 



2.33 2.7 7/2+ 



0.55(2) 3.05(3) 



0.40(4) 3.01 0.48 9/2+ 



3.63(9) 6.13(9) 0.26±g; 4 6 0.8(1) 6.25 0.94 5/2+ 
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t — i — i — i — | — i — i — i — i — | — i — i — i — i — | — i — r 

l9 C(p,p') E p =10 MeV; 




20 40 60 



C .m. ( de g) 

Fig. 2. Differential cross sections leading to the 1.46 MeV state in 19 C are compared 
to DWBA predictions obtained by assuming different initial-state configurations. 
.7^=5/2^ is assumed for the excited state. 
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10 1 



r( g .s.) 




20 40 



9 c.m. ( de g) 

Fig. 3. Differential cross sections leading to the 2.20 and 3.05 MeV states in 17 C 
are compared to DWBA predictions obtained by assuming different initial-state 
configurations. J n values assumed for the excited states are 7/2+ and 9/2+ for the 
2.20 and 3.05 MeV states, respectively. 
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